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MICROWAVE EMISSION FROM SPINNING DUST IN CIRCUMSTELLAR DISKS. 

Roman R. Rafikov 1,2 

Draft version February 5, 2008 

ABSTRACT 

In the high density environments of circumstcllar disks dust grains arc expected to grow to large 
sizes by coagulation. Somewhat unexpectedly, recent near-IR observations of PAH features from disks 
around Herbig Ae/Be stars demonstrate that substantial amount of dust mass in these disks (up to 
several tens of per cent of the total carbon content) can be locked up in particles with sizes ranging 
from several to tens of nanometers. We investigate the possibility of detecting the electric dipole 
emission produced by these nanoparticles as they spin at thermal rates (tens of GHz) in cold gas. We 
show that such emission peaks in the microwave range and dominates over the thermal disk emission 
at v < 50 GHz typically by a factor of several if > 5% of the total carbon abundance is locked up 
in nanoparticles. We test the sensitivity of this prediction to various stellar and disk parameters 
and show that if the potential contamination of the spinning dust component by the free- free and/or 
synchrotron emission can be removed, then the best chances of detecting this emission would be in disks 
with small opacity, having SEDs with steep sub-mm slopes (which minimizes thermal disk emission 
at GHz frequencies). Detection of the spinning dust emission would provide important evidence for 
the existence, properties, and origin of the population of small dust particles in protoplanetary disks, 
with possible ramifications for planet formation. 

Subject headings: accretion, accretion disks — circumstellar matter — stars: pre-main-sequence 



1. INTRODUCTION 

Recent IR and sub-mm observations of disks around 
young stars strongly support the idea that dust grains 
in these disks are significantly evolved compared to their 
ISM counterparts. There is mounting evidence that in 
many young stellar systems dust grains have sizes reach- 
ing cm (Testi et al. 2003; Acke et al. 2004; Natta et al. 

2004) , in contrast to sub- micron sizes of the ISM dust 
particles, which is a natural consequence of grain growth 
by coagulation. This process is often accompanied by 
the change of the material properties of dust grains indi- 
cated by their increased chrystallinity (van Boekel et al. 

2005) and by settling of dust towards the disk midplanc 
inferred from the change of the shape of the spectral en- 
ergy distribution (SED) of disks (Acke et al. 2004). Ob- 
servational evidence for dust growth exists not only in 
protoplanetary disks around T Tauri stars (stellar mass 
M* < 2 Mo), but also in their analogs around higher 
mass (2 M© < M„ < 10 M ) Herbig Ae/Be stars (van 
Kerckhoven, Tielens, & Waelkens 2002; Sloan et al. 2005; 
Schiitz, Meeus, & Sterzik 2005; Habart et al. 2005), as 
well as in disks around substellar objects (Apai et al. 
2005). 

At the same time, recent observations of the polycyclic 
aromatic hydrocarbon (PAH) spectral features in disks 
of Herbig Ae/Be stars (Meeus et al. 2001; Acke & van 
den Ancker 2004) provide evidence for the presence of 
a significant population of very small dust grains (sizes 
a ~ 3 — 100 A) in these disks. Previously, the exis- 
tence of a significant amount of carbonaceous nanopar- 
ticles in the ISM has been proposed to explain the IRAS 
observations of "unidentified infrared" emission features 
and strong mid-infrared emission component resulting 
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from the starlight reprocessing by the ultrasmall grains 
(Boulanger & Perault 1988). The proposed fraction of 
the ISM carbon mass locked up in very small grains 
considerably exceeds that implied by the extension of 
the conventional MRN dust size distribution (Mathis, 
Rumpl, & Nordsieck 1977) to very small sizes, below 50 
A. This suggests the existence of a separate population of 
very small carbonaceous grains which is distinct from the 
MRN distribution and contains ~ 10% of C in the ISM 
(Leger & Puget 1984; Draine & Anderson 1985). Thus, 
it is not surprising that nanoscale dust particles should 
exist at some level in disks 3 around Herbig Ae/Be stars 
since the latter were originally part of the ISM. However, 
the detailed modeling of the IR emission features in these 
systems suggests that PAHs may contain as much as sev- 
eral tens of per cent of total C mass in disks (Habart, 
Natta & Kriigel 2004). Thus, it appear that small dust 
population in these disks does not only withstand co- 
agulation, but on the contrary may have gained some 
additional mass, presumably due to the fragmentation 
of larger particles. This suggests very interesting size- 
dependent evolution of the dust population in Herbig 
Ae/Be disks. In the case of disks around T Tauri stars 
and brown dwarfs it is difficult to make definite state- 
ments about the presence of small dust grains because of 
the lack of strong stellar UV fluxes needed for exciting 
PAH molecules (although see Giirtler et al. 1999). 

In this paper we describe a new observational chan- 
nel for probing the presence and properties of very small 
dust grains in disks around young stellar objects. In 
the outer parts of these disks, ~ 10 2 AU from the cen- 
tral star, nanometer-size dust particles spin at thermal 

3 Possible origin of observed PAH emission in the roughly spher- 
ical envelopes around stars rather than in disks is unlikely given 
the correlation between the strength of the PAH emission and the 
disk geometry inferred from the SED shape (Habart et al. 2004). 
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rates of several tens of GHz. Because of their intrin- 
sic dipole moments, these "macromolecules" emit electric 
dipole radiation in the microwave band where the ther- 
mal Rayleigh-Jeans emission of the circumstellar disk is 
rather weak. This makes it possible to disentangle the 
spinning dust contribution to the disk emission and study 
the properties of this dust component. We describe the 
physics of the spinning dust emission and analyze its ob- 
servational signatures in &\ Applications of our results 
are discussed in 

2. EMISSION FROM SPINNING DUST GRAINS. 

Ferrara & Dettmar (1994) proposed that charged ISM 
dust grains spinning at thermal rates can produce ob- 
servable microwave emission via electric dipole radiation. 
This mechanism has been further investigated in detail 
by Draine and Lazarian (1998b; hereafter DL98) who 
considered, among other things, intrinsic dipole moments 
of grains, various processes governing the nonthermal 
spin rates of grains, and a size distribution of emitting 
nanoscale particles consistent with previous detections 
of PAH emission from the ISM. In this paper we will 
largely follow the approach adopted in these pioneering 
investigations. 

2.1. Circumstellar disk model. 

We describe the structure of the circumstellar disk us- 
ing conventional two-temperature flared disk model of 
Chiang & Goldreich (1997) which states that the vertical 
extent of the disk can be split into two well-defined re- 
gions: upper exposed layer in which dust absorbs most of 
the incoming starlight and reradiates it in the infrared, 
and lower shielded region which is hidden from direct 
starlight and is warmed only by the reprocessed radiation 
of the outer layer. For our purposes both layers can be 
considered isothermal with different temperatures. Tem- 
perature of the shielded layer T s h is lower than that of 
the overlaying exposed layer by a factor of several. 

The amount of mass contained in the exposed layer 
is very small: fraction of the total disk surface den- 
sity £ corresponding to this layer is ~ 7/[kp(T*)£], 
where 7 < 0.1 is the disk flaring angle and np(T+) is 
the Planck opacity at the stellar temperature. Since 
it is likely that even in the distant parts of the disk 
kp(T*)£ 3> 1, shielded region must contain most of 
the disk mass. Thus, if small particles are well mixed 
throughout the disk (which is a reasonable assumption 
for very small grains having settling times longer than 
the disk lifetime) most of the small dust particles must 
be shielded from the direct stellar illumination by the 
outer exposed layer. This inference significantly simpli- 
fies our subsequent treatment. 

We assume a simple power-law profile for the disk sur- 
face density 



scales as a power law of v and we adopt in this study 
(following the notation of Beckwith et al. 1990) 



£(r) - £irr Q , 



(1) 



where r% is the distance from the star normalized by 1 
AU, £1 is the surface density at 1 AU, and a is a pa- 
rameter that we can vary. Midplane temperature profile 
T s h(r) needed for the calculation of the grain spin rates 
is computed in Appendix 1X1 according to the prescription 
of the flared disk model of Chiang & Goldreich (1997). 

At the low temperatures of interest for us opacity k v 
is due to dust grains. At long wavelengths k„ typically 
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where z/12 = ^/(10 12 Hz) is the radiation frequency nor- 
malized by 10 12 Hz and K12 is the opacity at 10 12 Hz. 
Power law index /3 can be directly measured from the 
disk SED in the submillimeter range and in protoplan- 
etary disks it varies from 0.3 to 1.5 (Kitamura et al. 
2002), while ISM emission is characterized by (3 ps 1.7 
(Finkbeiner, Davis, & Schlegel 1999). Small values of (3 
in protoplanetary disks are interpreted as the evidence 
for grain growth (Natta et al. 2004; Draine 2005). 

Normalization of opacity at a given frequency K12 is 
rather poorly constrained. Beckwith et al. (1990) advo- 
cate using K12 =0.1 cm 2 g _1 while Kramer et al. (1998) 
find k(1.2 mm) ps 0.004 cm 2 g _1 which translates into 
K\i ps 0.016 cm 2 g _1 for (3—1. Draine (2005) modeled 
absorbing properties of evolved dust populations com- 
posed of different materials with MRN size distributions 
extending to maximum size a max ~ 1 — 10 cm and found 
ps 1 - 1.5 and k 12 ~ 0.01 - 0.1 cm 2 g" 1 (for the dust- 
to-gas ratio 10 -2 ). Thus, there is at least an order of 
magnitude spread in the possible values of K\2- 

2.2. Spinning dust emissivity. 

Midplane regions of circumstellar disks have very low 
degree of ionization because of intrinsic low temperature 
and strong shielding of the ionizing stellar radiation and 
cosmic ray fluxes by the overlying disk layers. As a result, 
dust grains in the shielded layer are predominantly neu- 
tral and their dipole moments are intrinsic and not due 
to an asymmetric charge distribution. Following DL98 
we assume very small dust grains to be composed of ran- 
domly oriented chemical substructures so that the total 
intrinsic dipole moment of the grain scales with the num- 
ber of atoms in the grain N as 



N^do. 



(3) 



Based on the available laboratory data on the dipole 
moments of the PAH-like particles DL98 have chosen 
do = 0.4 debye in their study, and we adopt this esti- 
mate as well. 

The power emitted by a grain of radius a spinning at 
frequency u> is given by 



4 d 2 {a)Lj 4 
9 c3 



167r d^uj pa 



27 



fid 



(4) 



where p is the bulk density of grain material and fid is the 
mean mass of grain atom (fid ps 9.25 amu for C:H=3:1 
carbonaceous material typical for PAHs). In deriving 
equation Hi we assumed that grains are spherical [so that 
d 2 (a) = do\4:-K pa 3 ) / (3fid)} and that their dipole moments 
are randomly oriented with respect to their rotational 
axes. 

To calculate the spectrum of the dipole emission of 
spinning grains we also need to know the size distribution 
of nanoscale dust particles. Contrary to the expectation 
of small dust depletion in circumstellar disks, recent ob- 
servations of PAH features in Herbig Ae/Be disks suggest 

4 DL98 assumed that in dense interstellar clouds PAHs are de- 
pleted by a factor of 5. 
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that very small grains may actually be more abundant 
in these systems than in the ISM (Habart et al. 2004). 
By analogy with DL98, we use the following functional 
form for the grain size distribution 



1 dn 
nn da 



= A 



MRN" 



-3.5 



B 



■ exp 



ln(a/a ) 



(5) 



where A M kn = 6.9 x 10~ 26 cm 2 5 (Draine & Lee 1984) 
and B are constants, ao is the typical size of very small 
carbonaceous grains, and a is the width of their log- 
normal distribution. The first term is a conventional 
MRN contribution while the second accounts for very 
small grains. Coefficient B is chosen in such a way that 
the fraction of C locked up in very small grains (with re- 
spect to the total abundance C/H=4 x 10~ 4 ) is fc- Our 
baseline model uses fc = 0.05 (close to the ISM value), 
ao = 3 A and a — 0.5, this distribution extends down to 
a minimum size a m in = 3.5 A. 

2.3. Grain spin rates. 

Rotation rates of small dust grains are in general gov- 
erned by a combination of physical processes: collisions 
of grains with neutral molecules, Coulomb interactions 
with charged species, infrared and electric dipole emis- 
sion, etc. (DL98). Since here we are mainly interested in 
obtaining reasonable estimate of the importance of the 
spinning dust emission in disks, we assume that grain 
rotation is determined by thermal equilibrium. 

This assumption is reasonable in the midplane region 
of the circumstellar disk because of the high gas density 
and very low ionization fraction there. Consequently, in- 
teractions of predominantly neutral dust grains with ex- 
tremely rare charged particles (plasma drag/excitation 
and rotational excitation by ion collisions, which are of- 
ten the most important determinants of grain spin rates 
in the ISM, see the case of molecular cloud environment 
in DL98) are completely negligible in this part of disk. 
As a result, grain- neutral collisions bring particle spin 
rates in thermal equilibrium with the surrounding gas 5 . 
This implies that grain spin rates have Boltzmann dis- 
tribution 

3/2 9 

■ exp 



with the rms rotation rate [uj 2 ) 1 / 2 given by 



(uj 2 ) 1 ' 2 = 



1/2 



3.5 x 10 10 T 2 1/2 a% 2 



(6) 



(7) 



where T-i is the gas temperature normalized to 10 2 K, 
a_7 is the grain radius in units of 10 -7 cm, and we take 
p = 2 g cm~ 3 . For simplicity, we assume grain shapes 
to be close to spherical implying that their moment of 
inertia is / = (87r/15)pa 5 , although very small particles, 
like PAHs, are likely to be sheetlike (see DL98 for more 
refined treatment). 

Within the upper exposed layer grains may spin at 
significantly non-thermal rates (corresponding to tem- 
peratures higher than the midplane temperature) but 
the amount of mass contained in this layer is very small 
meaning that its contribution to the total spinning dust 
emission can be neglected in the first approximation. 

5 Infrared emission by grains is unlikely to strongly affect their 
spin rates because of low midplane temperatures. 
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Fig. 1. — Spectrum of the microwave disk emission for systems 
described by Models 1 (a), 2 (6), and 3 (c) placed 100 pc away 
from us. Rayleigh- Jeans tail of the thermal disk emission {dashed 
line), contribution of the spinning dust emission {dotted line), and 
their sum {solid line) are displayed. Thermal emission is calculated 



for the dust emissivity index j3 = 1 and behaves as F* 



(see 



text for other parameters) . At radio wavelengths the spinning dust 
emission exceeds the thermal disk emission by F£ d / ' F* h ^ 3 (a), 
w 6 (6), « 4(c). 

2.4. Spectrum of the spinning dust emission. 

Combining equations ©, and JBJ we find the total 
spectrum of the spinning dust emission produced by the 
disk: 



Ff, d = 2tt 



1 dn 
da — 

nn da 



Rout 



■ I ^ {r)rdr P{a,w) x2nf u (u) 



1/2 , ,6 

= 8^ 2 f^ 1 



3tt 



nine 



/ —d (a) da 

J n H da 



T,(r)rdr 



exp 



(8) 



where Ri n and R ou t are the inner and outer radii of the 
disk, (u> 2 ) is given by equation Q as a function of T(r), 
and the maximum grain size a max is irrelevant for this 
calculation since F* d is dominated by small grains. 

For a simple case of a power law temperature distri- 
bution T(r) = To(r/i?*) -7 and a dust population with a 
single characteristic size a one finds that 

Vout <v<v in , (9) 



F sd K 1/ 3+2(a-2)/ 7 



where 



^out, in 



( 15 k B T 
\47r£ pa 5 



1/2 



ut,in 



-7/2 



(10) 
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are the frequencies determining the extent of the power 
law segment of F£ d . Outer regions of irradiated disks can 
often be well characterized by 7 ~ 0.5 which yields slope 
of the power law section of F£ d quite different from 2 + (3 
characteristic for the thermal disk emission 6 

Rout 

F? = 2tc J B v {T(r))ip sh (T{r))rdr (11) 

Rin 

in the optically thin Rayleigh- Jeans regime (here B y is 
a Planck function and tjj s h is defined in Appendix 0) : 
power law index of F£ d is —1 for a = 1 and 1 for a = 
3/2). As a result, one can hope to detect the spinning 
dust emission only at long wavelengths where F£ d may 
dominate over F* h . 

In Figure ^ we plot spectra of the spinning dust emis- 
sion computed according to the prescription (JHJ) for cir- 
cumstellar disks with different central objects at different 
evolutionary stages - massive Herbig Ae/Be star, solar- 
type T Tauri star, and a brown dwarf. Adop ted param- 
eters of these systems are listed in Table IA1I and are the 
same 7 as in Dullemond & Natta (2003). In all models 
we use a = 1, (3 = 1, n 12 = 0.05 cm 2 g" 1 and adopt a 
baseline model of the grain size distribution © with pa- 
rameters listed in 2.21 Disks are assumed to be located 
100 pc away from us. 

One can see that all three models exhibit significant ex- 
cesses due to the spinning dust emission around 30 GHz. 
This component exceeds the thermal power law contri- 
bution by a factor ranging from 3 to 6. Excess appears 
first at frequencies of about 50 GHz - typical spin fre- 
quency of dust grains at low temperatures, see equation 
J7J - and extends all the way to lower frequencies. 

2.5. Sensitivity to stellar and disk parameters. 

The appearance of the SED including the spinning dust 
component depends on a large number of stellar, dust, 
and disk parameters: M+, i?*, Si, a, f3, K12, etc. To test 
sensitivity to different parameters we have varied them 
one at a time to see how spectrum evolves with respect 
to a fixed fiducial model of the disk+star system. As a 
fiducial stellar model we have adopted T Tauri model 2 
but with M* = 1 M* and ai n — 0.3 AU for the inner 
edge of the disk. Fiducial model also uses a — 1, /3 = 1, 
K12 = 0.05 cm 2 g _1 . Parameters of the fiducial dust 
model a re sp ecified in i)2 2.21 

Figure lA"2l demonstrates how the total microwave spec- 
trum of the disk F£ d + F* h changes due to the param- 
eter variation. Figure IA3I shows the same but for the 
ratio of fluxes F^ d /F^ h . These plots demonstrate that 
the microwave spectrum is insensitive to either stellar 
mass M* or Ri n . The former appears only in the cal- 
culation of the radial temperature profile (see Appendix 
A) which depends on M+ only weakly. The inner radius 
of the disk does not play any role since at cm wave- 
lengths both the spinning dust emission and the thermal 
emission of the disk are dominated by the outer region, 

6 This expression accounts for the fact that at long wavelengths 
F^ h is dominated by the cool shielded region and the contribution 
of the warm exposed layer is small. 

7 Except the surface density normalization in the T Tauri model 
which in our case is comparable to the Minimum Mass Solar Nebula 
value. 



T ~ Rout- Variation of R ou t changes mass of the emit- 
ting gas which mainly affects the overall flux level (Fig- 
urc IA2b 'l but also changes F£ d /F* h somewhat (Figure 
IA3b ). The same is true for the spectral sensitivity to 
Si and a - variation of both parameters mainly affects 
the emitting dust mass, so the characteristic pattern of 
spectral evolu tion is simi lar to that caused by changing 
Rout (Figures lA2b .f and IA3b .f). Spectral sensitivity to 
T* and i?* is similar as both parameters determine the 
stellar luminosity. Changing the parameters of the small 
dust size distribution - do, cr, and fc - affects both shape 
and normalization of the spinning dust spectrum leaving 
the thermal disk emission (dominated by large grains) 
unaffected. 

Variations of K12 do not affect the spinning dust emis- 
sion but change the op tical depth of the disk which then 
affects see Figure lA"2l . As a result, for low enough 
opacity (e.g. K12 = 0.01 cm 2 g _1 ) spinning dust emission 
can dominate over the thermal disk emission by more 
tha n an order of magnitude at v < 30 GHz (see Figure 
IA3I ). Similar effect is produced by the variation of the 
dust emissivity index (3 - it changes the thermal flux of 
the dis k by orde rs o f magnitude since FfJ 1 cx v 2+l3 , see 
Figures lA2k andlAlk. 

Total power produced by the disk determines whether 
it can be detected in a flux limited observation. But 
even if such a detection is made, one can disentangle the 
spinning dust component from the thermal disk emission 
only if F£ d /F 1 *' 1 > 1. Examination of Figure ^3] demon- 
strates that this is most easily done for small K12 and/or 
large /?. Coagulation of dust particles in dense circum- 
stellar disks environments generally tends to decrease (3 
(sometimes below 1), which would hide the spinning dust 
emission bump in the thermal power law compo nent of 
the spectrum (see the curve for (3 = 0.5 in Figure lA"2k ) . 
Thus, from the observational point of view the most fa- 
vorable chances of detecting the spinning dust emission 
would be in systems having steep sub-mm slopes of the 
SED {(3 > 1). 

3. DISCUSSION. 

There is strong evidence for the presence of the spin- 
ning dust emission in the data obtained by the CMB ex- 
periments (Draine & Lazarian 1998a; Finkbeiner 2004; 
although see Bennet et al. 2004 for the opposite con- 
clusion) and by the recent Green Bank Galactic Plane 
Survey (Finkbeiner, Langston, & Minter 2004). It ap- 
pears that the microwave emission of the spinning dust 
in the ISM is an established physical phenomenon with 
substantial observational support. For this reason there 
is little doubt that it should also be present in the 
circumstellar disk environments and the two remain- 
ing issues would then be of quantitative nature: what 
is the strength of this emission and whether it could 
be observed. Results presented in the previous section 
clearly demonstrate that the spinning dust contribution 
to the microwave emission is substantial and in many 
cases significantly exceeds the thermal dust emission at 
v < 30 - 50 GHz. 

The answer to the second question strongly depends 
on the presence of other spectral contributions in the mi- 
crowave band. These include free-free emission from the 
gas heated by the central star or by the nearby massive 
stars and synchrotron emission from jets and outflows 
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which accompany large fraction of the young stellar ob- 
jects. Spinning dust emission can be discriminated from 
these components based on its intrinsic bell-like spectral 
shape since both free-free and synchrotron are character- 
ized by power-law spectra. Power-law spectral tails at cm 
wavelengths extending sometimes for more than a decade 
in frequency have indeed been observed in some objects 
(Zapata et al. 2005; Eisner & Carpenter 2006) suggesting 
the contamination by free-free and synchrotron emission. 
Given broad enough frequency coverage in the radio one 
can detect these power law components at lower frequen- 
cies, extrapolate them to 10 — 50 GHz and subtract from 
the total flux. Spinning dust emission will be present in 
the remaining flux provided that it dominates over such 
power law contaminants (its detection would be facili- 
tated by its intrinsic spectral shape). 

Another way of deducing possible free-free contami- 
nation is to look for Ha emission from the same object 
since both components originate in the same hot gas. 
This would give a useful handle on the expected strength 
of the free-free flux at low frequencies and demonstrate 
a priori whether it can compete with the spinning dust 
emission in the microwave domain. Synchrotron con- 
tribution can be identified based on its polarization and 
the morphology of the synchrotron emitting outflows (for 
which reasonable quality imaging would be required). 

Provided that these difficulties can be successfully 
overcome, the potential detection (or nondctection) of 
the spinning dust emission can give us very interesting 
information on the small dust population in the circum- 
stellar disks. There is already strong evidence for the ex- 
istence of such population in disks around Herbig Ae/Be 
stars based on the detection of PAH emission features in 
the 6 — 12 [im band (van Kerckhoven et al. 2002; Shiitz 
et al. 2005; Sloan et al. 2005; Habart et al. 2005) excited 
by copious UV flux of the central star. The total mass of 
nanoparticles producing near-IR emission features is very 
small (10~ 4 - 10 -5 M e , see van Kerckhoven et al. 2002) 
but one has to remember that the favorable PAH exci- 
tation conditions exist only in the very thin outermost 
layer of the disk where the stellar UV flux gets fully ab- 
sorbed. As a result, emission in PAH bands shows only 
the very tip of the iceberg as it traces an insignificant 
fraction of the total population of small dust. Most of the 
nanoparticles must be hidden from stellar UV near the 
disk midplane (provided that they are homogeneously 
distributed) and constitute substantial fraction of the to- 
tal dust population (containing up to 30% of the total C 
abundance, see Habart et al. 2004). As the spinning dust 
emission is contributed by all small dust particles (disk 
is optically thin at microwave frequencies), its detection 
would give us a very useful probe of the total mass in 
nanoparticles. Combined with the model for the disk 
structure, the comparison of the near-IR PAH observa- 
tions and the spinning dust emission would provide a 
very good consistency check for the two probes of the 
small dust population. 

A nondetection of the spinning dust emission in Herbig 
Ac/Be systems may suggest that either (1) vertical mix- 
ing in disks around these stars is for some reason inhib- 
ited (which is unlikely given that the magnetorotational 
instability is expected to operate at least in the upper 
disk layers) or (2) PAH formation is catalyzed by stellar 
UV in the uppermost disk layer and they are rapidly de- 



stroyed in its interior. Both possibilities appear rather 
exotic. Note that because of the high stellar tempera- 
ture the contamination with free-free emission may be 
especially severe in Herbig Ae/Be systems, potentially 
allowing only the upper limits to be set on the flux due 
to the spinning dust. 

With only few exceptions (e.g. Giirtler et al. 1999) 
mid-IR PAH emission has not been detected in disks 
around T Tauri stars or brown dwarves. This fact is 
consistent with the lack of strong UV flux in this ob- 
jects, which is necessary for efficient excitation of the 
PAH features (although see Li & Draine 2002 and Smith, 
Clayton, & Valencic 2004 for a different opinion), and 
should not be immediately interpreted as the evidence 
for the absence of the nanoparticle dust component in 
these systems. Because of lack of the PAH emission, the 
microwave emission of spinning dust may be the only 8 
way in which the small dust could be probed in T Tauri 
and brown dwarf disks. 

Unlike the mid-IR PAH features, spinning dust emis- 
sion reveals the whole nanoparticle dust population of 
the disk, which is important in many respects. Small 
particles affect the ionization balance within the disk 
since they dominate the dust surface area and are very 
efficient at immobilization of free charges. This works 
against the operation of the magnetorotational insta- 
bility in protoplanetary disks (Fleming & Stone 2003) 
which, in turn, affects the sizes and positions of the 
so-called "dead zones" of quiescent non-turbulent ma- 
terial (Gammie 1996), with potential consequences for 
the planet formation (Matsumura & Pudritz 2005, 2006). 
Small grains also affect heating of gas in the upper rar- 
efied disk layers due to the higher yield of electrons from 
PAHs (Kamp & Dullemond 2004). 

The very existence of substantial amounts of nanopar- 
ticles is very interesting from the evolutionary point of 
view since dust grains are expected to grow as disks age. 
Recent discoveries of PAHs in Herbig Ae/Be disks in 
amounts exceeding those in the ISM (Habart et al. 2005) 
are at odds with this expectation but may be in line 
with the recent work of Dullemond & Dominik (2005) 
who found that dust fragmentation, in addition to co- 
agulation, has to be an important ingredient of the dust 
evolution. It may also be the case that the material prop- 
erties of small and large grains (e.g. sticking coefficient) 
are completely different leading to divergent evolutionary 
paths of the two extremes of the grain size distribution. 
Detection of the spinning dust emission from circumstel- 
lar disks may shed enough light on the properties of the 
small dust population to resolve this puzzle. 

Calculations presented in fJ21 w CTe intended to provide 
us with the rough idea for the importance of the spin- 
ning dust emission and as such they neglected a number 
of details. Among them are the realistic distribution of 
shapes of the nanoscale dust particles, emission of spin- 
ning dust in the exposed layer associated with the non- 
thermal grain rotation rates, possible ionization of small 
dust particles, and so on. In the future observational de- 
mands may warrant a more refined study including these 
details. 

8 The contribution of very small particles to opacity is degenerate 
with that of large grains. 
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4. SUMMARY. 

We have studied the electric dipole emission produced 
by spinning dust grains of very small size (several to tens 
nm) in circumstellar disks. At the temperatures char- 
acteristic for these environments spinning grains emit in 
the microwave range with the peak of the emission occur- 
ring at ~ 30 — 50 GHz. For typical parameters of disks 
around Herbig Ae/Be stars, T Tauri stars, and brown 
dwarves we find spinning dust emission to dominate over 
the thermal disk emission by a factor of at least several 
at v < 50 GHz (if nanoparticles contain > 5% of the 
total C abundance). We have studied the sensitivity of 
this emission component to various stellar, disk, and dust 
parameters and found it to be strongest for variations of 
the total dust opacity — both normalization of k„ and 



the slope of its frequency dependence. Provided that 
the contamination of the spinning dust emission by the 
free-free and/or synchrotron emission can be mitigated, 
the best chances of detecting this emission component 
should be in systems with steep sub-mm spectral indices 
minimizing the contribution of the thermal disk emis- 
sion. Detection (or non-detection) of the spinning dust 
emission will provide important information about the 
existence, properties, and origin of the nanoscale dust 
particles in circumstellar disks. 



RRR thankfully acknowledges the support of this work 
by the Canada Research Chairs program and Connaught 
Foundation. 



APPENDIX 

RADIAL TEMPERATURE PROFILE IN THE DISK MIDPLANE. 

Thermal balance of the shielded region is set by (S / 4)Tj (R+ / a) 2 ip ex = T^ipsh, where T sh and T ex are the character- 
istic temperatures of the shielded and exposed regions respectively and ip ex and ip s h are the fractions of the radiative 
flux with blackbody temperatures T ex and T s h absorbed by the disk. In the optically thick case ipex,sh = 1 while in the 
optically thin regime ip e x,sh = ^>Kp(T eXiS h). Flaring angle <5 = d(Hi/a)/dhia is determined by the height of the disk 
surface Hi w Ac s (T s fc)/fi, where c s is the gas sound speed and f2 is the angular frequency of the disk. Factor A ~ 1 is 
roughly constant throughout the disk and we set A = 3 in this study. Power law scaling of k„ implied by equation J3J) 
results in the following expression for the low temperature Planck opacity: 

K P (T) = K* — , k* = —_— — j — / — - (Al) 



7T 4 (10 12 Hz)^ 



- 1 



where (3 is the same as the power law index in equation © . Opacity dependence given by (|A1|) applies to the emission 
of both the disk interior and the exposed layer, but not to the stellar radiation, so that Kp(T*) is significantly different 
from In our case Kp(T\) is calculated for a population of small 0.1 /im dust grains (Draine & Lee 1984; Dullcmond 
& Natta 2003). Using (|AT) one finds that 

- Kp (T^-\ V(4+/3) /D \ 2/(4+0) 



T ex {r) = 



4k* 



J2* 



(A2) 



This information allows one to derive the temperature profile in the shielded layer in different parts of the disk. 
Inner disk is optically thick to radiation at both T ex and T sh , which results in 



T sh (r) — T+ ( 
shy j *\ UR ^ 



2/7 



r 

R~* 



-3/7 



In the intermediate region disk is optically thin to the radiation of its interior [ip s h — Y>Kp(T s h 
thick to the reprocessed emission of the exposed layer (ip ex = 1). As a result, 



T sh (r) =T* 



2/(7+20) , s 

where = S(i2*) is a value of E obtained by extrapolation of QJ to r 
optically thin to the radiation at both T ex and T s h, one finds 



2 + a + (3 A h+ 
4(7 + 2/3) k^R~* 



(2a-3)/(7+2/3) 



(A3) 

< 1] while still optically 
(A4) 



R+. Finally, in the outer disk, which is 



T sh (r) = T* 



(3 2 + 4/3 + 8 
4(4 + /3)(7 + 2/3)"i?* 



a^ 



2/(7+2/3) r 



Kp(T^) 



4k* 



2/3/[(4+/3)(7+2/3)] 



r 

R~* 



-(7/3+12)/[(4+/3)(7+2/3)] 



Transition between the inner and intermediate regions occurs at 

7/(7q+3/3) I A h* 



ai ps i?* (k*E*) 
and between the intermediate and outer regions at 

a 2 *R+( K ^f +Wl2f3+a(4+m 



14 R* 



2/3/(7a+3/3) 



| 3/[2/3+ct(4+/3)] 



(A5) 



(A6) 



(A7) 



In a particular case of a — 3/2, (3 = 1 our results a gree with Chiang & Goldreich (1997). The profile of T s h(r) used 
in [2| is found by extrapolation between expressions ljA3(l - l|A5(l . 
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Fig. A2. — Evolution of the total disk spectrum (both thermal and spinning dust contributions included) caused by the variation of 
basic model parameters. Each window displays the results obtained by variation of a single parameter around the fiducial model described 
in the text. The parameter being varied and its values corresponding to different curves are indicated in each window. Solid curve always 
represents the spectrum of the fiducial model. 



9 




Fig. A3. — Evolution of the ratio F£ d /F* h representing the strength of the spinning dust contribution relative to the the rmal disk 
emission due to the variation of the model parameters with respect to the fiducial model. Notation is the same as in Figure lA*2l 



